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The formation of Bose condensates far from equilibrium can play an important role in our under- 
standing of collision experiments of heavy nuclei or for the evolution of the early universe. In the 
relativistic quantum world particle number changing processes can counteract Bose condensation, 
and there is a considerable debate about the relevance of this phenomenon in this context. We show 
that the involved question of Bose condensation from initial over-population can be answered for 
the example of scalar field theories. Condensate formation occurs as a consequence of an inverse 
particle cascade with a universal power-law spectrum. This particle transport towards low momenta 
is part of a dual cascade, in which energy is also transfered by weak wave turbulence towards higher 
momenta. To highlight the importance of number changing processes for the subsequent decay of 
the condensate, we also compare to non-relativistic theories with exact number conservation. We 
discuss the relevance of these results for nonabelian gauge theories. 
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Seventy years after the prediction of Bose and Einstein 
for a new state of matter of a dilute gas of weakly cou- 
pled bosons the first gaseous condensate was produced 
using ultracold rubidium atoms [l]. By now a very de- 
tailed understanding exists, both experimentally as well 
as theoretically, about how a condensate emerges by de- 
creasing the temperature going through a sequence of 
close-to-equilibrium states. Much less is known about 
Bose condensation far from equilibrium. Most impor- 
tantly, this concerns open questions about condensation 
dynamics in relativistic quantum field theories includ- 
ing particle creation and annihilation processes. Particle 
number changing processes might change the situation 
dramatically as compared to the number conserving non- 
relativistic counterparts. This can play a crucial role in 
our understanding of far-from-equilibrium stages of the 
early universe [2[ as well as for related questions in rela- 
tivistic collision experiments of heavy nuclei 0] . 

Recently, it has been argued that Bose condensation 
in relativistic quantum field theories occurs for a generic 
class of far-from-equilibrium situations. For the exam- 
ple of a plasma of quarks and gluons in the early stages 
of a heavy-ion collision at sufficiently high energies, the 
gluon density may be parametrically large compared to 
the thermal equilibrium value As a consequence of 
the initial over-population, where the gluon occupation 
number at a large characteristic momentum scale Q s is 
parametrically of the order of the inverse gauge coupling 
l/a s (Q s ) 3> 1) condensation is suggested to occur. Sim- 
ilar situations arise also in the context of scalar infla- 
ton dynamics in the early universe following the violent 
process of preheating Q. Here a nonequilibrium insta- 
bility, such as parametric resonance in chaotic inflation- 
ary scenarios or a tachyonic instability, leads to over- 
population of typical modes with occupancies inversely 
proportional to a small quartic self-coupling A < 1. De- 



spite the weak coupling, the parametrically large occu- 
pancies lead to strong correlations such that the problem 
is non-perturbative. 

In this work we show that Bose condensation far- 
from-equilibrium generically occurs from initial over- 
population as a consequence of a particle cascade towards 
low momenta. We consider relativistic scalar field theo- 
ries where the occupation number at a large character- 
istic momentum scale, n(Q s ), starts out proportional to 
the inverse quartic self-coupling, 1/A> 1. Remarkably, 
as a consequence of this non-perturbative occupancy the 
condensation dynamics becomes universal. For a large 
class of models, independently of the value of the (small) 
coupling or mass or (sufficiently large) cutoff scale, the 
same phenomenon is observed. A dual cascade develops: 
for momenta p^Q s a power-law distribution 

n(p>Q s )-^f) (1) 

signals the direct energy cascade towards higher mo- 
menta known from the theory of weak wave turbulence 
in d spatial dimensions [a, @] . This energy cascade goes 
along with an inverse particle cascade towards the in- 
frared with the low-momentum distribution 

n(p<Q s )~(^y (2) 

known from strong turbulence applications to early- 
universe inflaton dynamics [7|, |8J . Corresponding power- 
law regimes with directed fluxes have also been found in 
non-relativistic systems of cold atoms 0, [T3| • 

We find that the transport of particles from the in- 
frared cascade leads to condensation of the zero momen- 
tum mode of the anti-commutator expectation value of 
two Heisenberg field operators <p{t,x), 
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(2PI) resummed large- N expansion to next-to-leading or- 
der for iV-component scalar field theories @, EH ■ It es- 
sentially sums elastic and inelastic processes to infinite 
order as a geometric series, which can be encoded into 
an 'effective coupling' X cfi (p) = A/|l + Ur(p)\ 2 . Here 
Hr(p) = ^Ir{p° = p,p) denotes the self-consistently 
dressed (2PI) one-loop retarded self-energy with on-shell 
frequency p° — p [T3|- Since Hr(p > Q s ) < 1 one finds 



(5) 



such that standard kinetic descriptions are recovered at 
large momenta. In contrast, using the results of Ref. 0], 
for small momenta Hr(p) 3> 1 such that the effective 
coupling encodes the remarkable infrared scaling 



FIG. 1: Schematic picture of Bose condensation far-from- 
equilibrium from the dual cascade. Also indicated are the 
ranges of validity of kinetic theory, classical-statistical field 
theory and the 'vertex' resummed kinetic theory based on 
the two-particle-irreducible (2PI) large- N expansion to NLO. 



In spatial Fourier space, F(t 7 t';p) = J d d p/(2n) d exp(ip ■ 
x) F(t,t'-x) with p = \p\ can be decomposed into a 
(time-dependent) condensate fraction ~ 0q (t) and a non- 
condensate fraction involving the distribution function 
n p (t) = n(t,p) as well as the dispersion cu p (t), 

F(t,t;p) = -1^ ("*(*) + 5) +(2n) d 6^(p)cbl(t)- 

The situation for Bose condensation far from equilib- 
rium is summarized schematically in Fig. [TJ Displayed 
are the characteristic regimes of the distribution as a 
function of momentum on a double logarithmic scale. At 
large momenta p 3> Q s , where occupancies n{p) < 1, 
a quantum or dissipative regime preempts a power-law 
behavior. At lower but still large momenta p > Q s , 
where parametrically the occupation number is in the 
range 1/A> n(p) ^> 1, the energy transport occurs via 
the direct cascade towards higher momenta. In particu- 
lar, the power-law exponent d — 3/2 can be understood 
as a consequence of the presence of a condensate leading 
to an effective cubic interaction In the low momen- 
tum regime p < Q s occupancies grow non-perturbatively 
large, i.e. parametrically n(p) ~ 1/A. Here the inverse 
particle cascade leads to condensation at p = 0, which is 
indicated schematically on this logarithmic plot. 

We emphasize that condensation from initial over- 
population is not described by conventional kinetic the- 
ory including a finite number of elastic or inelastic pro- 
cesses. Instead, an infinite series of processes, all be- 
ing of order one, has to be summed to get the leading 
contribution. Fortunately, there exists a non-trivial ex- 
ample where this can be computed directly in quantum 
field theory. It is based on the two-particle irreducible 
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Taking into account also the power-law behavior ([T} and 
([2|) of the distribution function, scaling analysis shows 
that elastic processes dominate the regime p < Q s - The 
non-condensate fraction of the distribution follows then 
the effective number conserving 'kinetic equation' 



^ ~ A / dfi 2o2 [A 



eff , \eff , xeff ] 
p+l ' A p-q ' A p-r\ 



x [(n p + ni)n q n r - n p ni(n q + n r )] , (7) 

where J ; r df22<-s-2 denotes the relativistically invariant 
measure with energy-momentum conservation for two-to- 
two scattering 12j . The infrared particle cascade ([2]) is a 
stationary solution, dn p /dt = 0, of this equation I n 
the regime p > Q s , where ([5]) holds, the energy cascade 
driven by particles scattering off the condensate frac- 
tion takes over. The corresponding kinetic description 
in this momentum regime can be formally obtained from 
([7]) with the replacement 



(2n) d S^(p)cPl 



(8) 



The energy cascade fij) is a stationary solution of the 
resulting equation [f| 0]. 

This scenario of far-from-equilibrium Bose condensa- 
tion can be rigorously tested using numerical simulations 
on a space-time lattice. Here we exploit the fact that 
classical-statistical field theory descriptions 13j have an 
overlapping range of validity with the above effective ki- 
netic theory 14| or the underlying quantum field the- 
ory Q. As a consequence, the dynamics of condensa- 
tion and the relevant cascading regimes with occupan- 
cies n(j>) 3> 1 can be accurately described by numerically 
solving the classical field equations of motion and Monte 
Carlo sampling of initial conditions. 

In Fig. [5] we show the result of simulations for a 
relativistic N — 4 component scalar field theory with 
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FIG. 2: Starting from an over-populated initial distribution, 
a dual cascade develops in the relativistic theory. 



A/(4!iV) J2 a =i('l )a< i )a ) 2 interaction in d — 3. Initial con- 
ditions are sampled to generate the non-zero classical av- 
erages (#)cl(* = 0,P) = Hp/"® and {4>4>) c \{t = 0,p) = 
ripUjp with initial dispersion u> p = \J p 2 + m 2 . Starting 
from a rescaled distribution n p X/6N = AQ(Q S — p) 
with amplitudes in the range A ~ 1 — 10 to obtain 
over-population up to the characteristic momentum scale 
Q s , we find that the system develops the dual cascade 
with the expected exponents during a long time interval 
tQ s ~ 100 — 10000. In particular, we find that the effec- 
tive mass scale m 2 /Q 2 <C 1 decreases with time while the 
energy and particle cascades operate. Our main result is 
given in Fig. [31 which shows the zero mode of the corre- 
lator F(t,p = 0)/V = (cj>4>) c i(t,p = 0)/V as a function 
of time for different volumes V. For finite V the Dirac 
(5-function in (j4} at zero spatial momentum is replaced 
by (27r) d j( d )(0) -4 V. Fig. |3] shows that initially this 
correlator is proportional to 1/V since the condensate 
fraction is zero. The subsequent emergence of a volume- 
independent contribution signals condensate formation. 

Particle number changing processes, which are present 
in the relativistic theory, finally lead to the decay of the 
condensate as can be seen from Fig. [3] Correspondingly, 
the energy cascade of the non-condensate fraction shown 
in Fig. [2] moves to the infrared at these times. In theories 
with exact number conservation a different behavior is 
expected. For comparison, we consider in the following as 
a well-known second example the non-relativistic Gross- 
Pitaevskii equation for scalar fields, which has been stud- 
ied extensively also in the context of ultracold atomic 
gases [l5(. The relevant coupling g — Ana/m in this case 
is characterized by the s-wave scattering length a and 
the mass m. Here we employ classical-statistical simula- 
tions of that equation to include fluctuations in complete 
analogy to the relativistic case above 0, [l(| • 

In Fig. 2] we show the results for the distribution n(p) 
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FIG. 3: Emergence of a volume-independent condensate frac- 
tion from runs with different volumes V = L s . 



in the non-relativistic theory starting from similar ini- 
tial conditions with over-population as above for d — 3. 
Again, the subsequent evolution builds up an inverse par- 
ticle cascade for p < Q s , which for the non-relativistic 
theory is characterized by the scaling exponent d + 2 [§] . 
Furthermore, we see that instead of the energy cascade 
at higher momenta rather quickly a classical thermal dis- 
tribution ~ 1/p 2 occurs. As a consequence of the particle 
transport towards lower momenta a condensate fraction 
appears, as is seen from Fig. [S] Total particle number 
/ d d p/(2ir) d n(p) is conserved and we observe no decay 
of the condensate due to number changing processes. 

We emphasize that the infrared particle cascade is 
an attractor solution, which is approached from a much 
wider class of initial conditions than the considered over- 
population scenario also in different dimensions 0, EH ■ 
Its physical nature may be associated to the formation of 
non-trivial topological configurations, which can be ob- 
served from the inhomogeneous field evolution before en- 
semble averaging [9l.[lol.fl6|. This provides an alternative 
viewpoint on the same physics as described by the above 
effective kinetic theory for homogeneous ensemble aver- 
ages. The latter also allows one to discuss cascades as 
non-equilibrium renormalization group fixed points [8| . 

Though the above scenario of Bose condensation far 
from equilibrium exhibits universal properties, which ap- 
ply to a large class of scalar field theories, its gen- 
eralization to theories with gauge bosons is far from 
trivial. Taking the relevant example of relativistic, 
nonabelian SU(N) gauge theories, for the perturbativc 
regime l/a s (p) ^> n s (p) 3> 1 it has been shown both 
from classical simulations and resummed perturbation 
theory that nonabelian gauge theories and scalar theories 
can exhibit the same class of weak wave turbulence expo- 
nents II 71 1 . This scaling behavior may be modified at later 
times [18| . In the kinetic regime one can study distribu- 
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FIG. 4: Similar to Fig. O however, for a non-relativistic the- 
ory described by classical-statistical Gross-Pitaevskii theory. 
The infrared particle cascade emerges with the predicted non- 
relativistic exponent d + 2 for the three-dimensional theory. 



FIG. 5: Similar as in Fig. [3l however, for the non-relativistic 
theory with conserved particle number. In contrast to the 
relativistic case, no decay of the condensate due to number 
changing processes is observed. 



tion functions which are derived from equal-time corre- 
lation functions, for instance, in Coulomb gauge. Us- 
ing classical-statistical SU{2) gauge theory simulations 
in d — 3 starting from over-populated initial conditions, 
we indeed find the same weak wave turbulence exponent 
3/2 as for the relativistic scalars presented above flil ]. 

However, it is unclear of how to interpret the most in- 
frared modes beyond the kinetic regime in this way and 
suitable gauge-invariant signatures associated to conden- 
sation are to be devised. A first non-perturbative study 
of coherent nonabelian gauge field dynamics is given in 
Ref. [20(. It has also been argued that the infrared cas- 
cading solutions observed for scalars Q can be carried 
over to nonabelian gauge theories [jlj . An infrared ex- 
tended kinetic description as given for the scalars above is 
complicated for gauge theories, since its straightforward 
generalization would involve the impractical resumma- 
tion of all planar diagrams. Classical-statistical simula- 
tions, which can also include the relevant physics of lon- 
gitudinal expansion in the context of heavy-ion collisions, 
should be the appropriate tool to clarify these questions 
for sufficiently high Q s . 
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